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CHAPTER I 


INTRODUCTION 


1.1 Introduction 

Penetrometers are widely used throughout the world by civil engineers for in-situ 
measurements of soil properties. A wide variety of penetrometers have been developed over 
the years, and the interpretation of results is therefore dependent on the type of penetrometer 
used and the testing procedure. In the classification of penetrometers, a distinction between 
Static and dynamic penetrometers is made. A static penetrometer uses a static force to 
penetrate an object into the soil, while dynamic penetrometers utilize an impact force for the 
penetration. The Standard Penetration Test (SPT) is one form of dynamic penetration test. 
The practice of driving a one-inch diameter open-end into the bottom of a bore hole to 
recover drive samples of soil was introduced in 1902 and make the beginning of dynamic 
penetration testing and sampling of soils. Between the late 1920’s and the early 1930’s, the 
attempt was made to standardize the test. The term "standardization" is a misnomer since 
many versions of the actual SPT system (hammer, rod and sampler) are being used 
throughout the country. 

Since the SPT blow count is inversely proportional to the energy delivered to the rod 
(Schmertmann and Palacios 1979), different methods and procedures have been suggested and 
used to calculate this energy so that the blow count could be adjusted to a standard energy. 


From wave mechanics theory, an accurate measurement of the energy transferred to the rod 
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can be obtained by integrating the product of the force and velocity measured at the top of 
the rod as a function of time during a hammer blow. In addition to the energy a force- 
displacement relationship is also available. 

Research in the pile driving area on pile capacity and pile-soil interaction has 
successfully utilized dynamic measurements of force and acceleration at the top of a pile 
during driving. This has been a major contribution to the advancement of the state-of-the-art 
in the pile dynamics area (Goble et al. 1975). From the records of force and acceleration, 
a prediction of the soil resistance forces along the pile was developed. Also, a calculation 
of the reaction and velocity at the pile toe was made from these records. 

Because of the obvious similarity between the Standard Penetration Test procedure 
and the mechanics of pile driving, attempts have been made to measure the dynamic event 
of a blow in the SPT, utilizing knowledge from the pile driving area. In this project, a 
measurement system similar to that used for piles was transferred to the hammer-rod system 
used in the SPT. This system was them tested in the laboratory. When good measurements 


were obtained the system was tested in the field during actual SPT operations. 


1.2 Scope and Objective 

The origin of the Standard Penetration Test, its applications to soil formations, the 
interpretation of data, and the limitations of the test are discussed in Chapter II along with 
a description of different procedures to measure the energy. A description of the calculation 
of force and displacement at the toe of the rod from the records of force and velocity at the 
top, is also given in this chapter. In Chapter III there is a theoretical discussion of wave 
propagation in SPT systems and a derivation of a means of calculating the wave form 


generated in the rod by a hammer blow is presented. Measuring systems, data acquisition 
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systems and processing procedures used for the tests in this project are described in Chapter 
IV. Evaluation of the results from laboratory and field tests are given in Chapter V, and the 
conclusions are summarized in the following chapter. 

This thesis describes the attempt to measure force and acceleration at a point in the 
rod during an SPT blow. No attempt has been made to study the sampler-soil interaction or 
to correlate SPT blow count to the engineering properties of soils. With an accurate 
measurement system of this kind, a better understanding and prediction of the above 


mentioned characteristics can be developed. 


CHAPTER II 


REVIEW OF PREVIOUS WORK 


2.1 Standard Penetration Test 

The Standard Penetration Test (SPT), defined by the American Society for Testing 
and Materials (ASTM) D-1586-67, is widely used in the United States for the investigation 
of subsurface soil properties. The test has developed from being an individually designed 
and locally used test into a nationally standardized and widely used test for soil exploration. 

The test consists of dropping a 140 lb hammer for a free fall of 30 inches to impact 
the top of a rod to which a split-spoon sampler is connected at the lower end (Fig. 2.1). The 
number of blows required to drive the sampler 12 inches into the ground is known as the 
standard penetration number, the N-value, or the blow count. The usual procedure is to 
drive the rod into the bottom of a borehole and to record the number of blows each 6 inch 
interval for the first 18 inches of driving. The blow count is taken as the sum of the last two 
6 inch increments. The reason for this is that the drilling operation is assumed to have 
disturbed the soil directly below the bottom of the boring, giving an erroneous blow count 


for about the first 6 inches. 
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Figure 2.1: A typical Split-Spoon Sampler by ASTM D1586-67. 


2.2 History of the SPT 

The Standard Penetration Test came into being as a result of the development of dry 
sample recovery techniques. Previously, subsurface investigations were performed primarily 
through the use of wash borings. A wash boring involves the circulation of a water and/or 
drilling mud mixture to remove the cuttings from the boring as the hole is advanced. In 
1902, Charles R. Gow introduced the first method of dry sample recovery (Fletcher 1965). 
He used a 110 Ib weight to drive a one inch outside diameter sampling pipe. After this 
method was used for a short time, it became apparent that the resistance to driving the 
sampler was influenced by the condition and properties (e.g., strength and density) of the 
soil, Thus, the term, “‘penetration resistance", was used to define the number of blows 
required to drive the sampler a given distance. 

In 1922, the Charles R. Gow Company merged with the Raymond Concrete Pile 
Company where L. Hart and G.A. Fletcher devised a split-spoon sampler (1927) of 2 inches 
diameter which is very similar to the current Standard Penetration Test (ASTM D-1586). In 
the late 1920’s through the efforts of Gordon Fletcher, Linton Hart and Harry Mohr, the first 
steps towards standardization were undertaken. The drive weight was changed from 110 Ib 
to 140 Ib and the drop height was set at 30 inches. After extensive field and laboratory 
work, the number of blows required to drive the sampler into the soil a distance of 12 inches 
was established as an official record of the test. After the initial work in the later 1920’s the 
test remained essentially unchanged until 1954 when James Parsons introduced a new method 
of recording the blow count. Instead of a penetration of 12 inches the sampler was driven 
18 inches into the soil and the blow count was noted for each increment of 6 inches. The two 
smallest 6 inch incremental blow counts were then added and recorded as the blow count. 


Since 1954, the only feature of the test to change significantly has been the method 
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of recording the blow count. The sampler is still driven 18 inches into the soil. The blow 
count for the first 6 inches of penetration is recorded but is considered a seating drive. The 
blows required to drive the sampler the final 12 inches are recorded as the blow count (N 
value) although the number of blows required for each six inch increment is normally 


recorded. 


2.3 Uses of the SPT 

The original purpose of the SPT was to obtain an approximate correlation between 
the blow count and expected behavior of the soil under load. However, as the popularity of 
the test increased, so did the uses in which it was employed. 

Researchers and engineers have made a great deal of effort to correlate the results of 
the SPT to the determinations of material properties and the development of design aids. The 
SPT is performed in granular and cohesive soil to obtain an indication of the in-situ soil 
properties. Terzaghi and Peck (1948) proposed a correlation between the N value, the 
relative density and bearing capacity of sand (Table 2.1); also a correlation between the N- 
value, the consistency and the unconfined compressive strength of cohesive soil (Table 2.2). 
Research on this subject was performed by Gibbs and Holtz (1957) and later by Marcuson 
and Bieganousky (1977) to estimate the in-situ relative density from the N-value. They also 
performed correlation studies between blow count, relative density of sand and the effective 
overburden pressure. See Figure 2.2. 

The blow count has also been shown to have a correlation with the shear strength 
parameters of soil (Fig. 2.3). De Mello (1971) and others found by statistical analysis of 
data that curves such as the ones in Figure 2.3 agree reasonably well with reality, provided 


the SPT blow count is not from very shallow depths. 


Table 2.1: Relative Density of Sand 


Blow Count Relative Allowable Bearing 


N density | Capacity gu(tons/ ft?) 


Less than 10 Loose Requires Compaction 


10 - 30 Medium 0.7 - 2.5 


30 - 50 Dense 2.5 - 4.5 


Over 50 Very Dense over 4.5 


Table 2.2: Consistency of Cohesive Soil 


Consistency | Unconfined Compressive 


Strength (tons/ ft?) 


very soit < 0.25 
soft 0.25 - 0.5 
medium 0.5 - 1.00 


stiff 1.00 - 2.00 


16 - 30 very stiff 2.00 - 4.00 


Over 30 hard over 4.00 
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Figure 2.2: Correlation Between SPT Blow Count and Relative Density of Cohesionless Soils 
as Given by Gibbs, Holtz and Bazarra. | 
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Figure 2.3: Estimation of @ (Angle of Internal Friction) from the SPT Blow Count. (After 
Schmertmann) 
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The test is also very useful in the design of pile foundations. It is used to simulate the 
pile during driving in an attempt to discover any problems that may occur when the actual 
driving takes place. Most pile foundations in granular soils are initially analyzed on the basis 
of SPT N-value. Meyerhof (1976) suggested that the ultimate point resistance in a 
homogeneous granular soil can be obtained from standard penetration resistance (Fig. 2.4); 
also he indicated that the average unit frictional resistance for driven high-displacement piles 
can be obtained from average standard penetration resistance values N (Fig. 2.5). In which 
N = N when N < 15 and N = 15+1/2(N-15) when N > 15. The test is very useful in 
determining in-situ soil coefficients commonly used in computer programs that use the one- 
dimensional wave equation to model a pile during driving. 

The Standard Penetration Test is also used for evaluating the liquefaction 
characteristics of a sand deposit. It was found that the most significant property of a soil 
affecting its liquefiability was the relative density, and as previously mentioned the SPT can 
give an indication of this. Schmertmann (1977) indicated such a relation, and Seed et al. 
(1977) found the correlation shown in Figure 2.6 where N-value was corrected to an 
overburden pressure of 1 ton/sq.ft, because of the effects of increasing overburden pressure 
on penetration resistance. 

In addition to the information obtained during the performance of the test the sample 
can yield some very useful data. Since the sample is highly disturbed, only identification 
tests are usually performed. From the sample the engineer can determine color, 
macrostructure, strata delineation, water content, Atterberg limits, specific gravity and 


texture. 
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Figure 2.4: Emperical Relation between Ultimate Point of Piles and Standard Penetration 
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Figure 2.5: Empirical Relation between Ultimate Skin Friction of Piles and Standard 
Penetration Resistance in Cohesionless Soil (1 tsf = 95.8 kn/m’) 
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Figure 2.6: Correlation Between Stress Ratio Causing Liquefaction in Field and Penetration 
Resistance of Sand. 


15 
2.4 Reliability of Blow Count of the SPT 

The reliability of the blow count depends on the ability to maintain consistent 
delivered energy in driving systems. Different delivered energies may result in significantly 
different blow counts in the same deposit at the same overburden pressure because the SPT 
blow count has been shown to be inversely proportional to the delivered energy 
(Schmertmann and Palacious 1979). An understanding of the factors which affect the 
penetration resistance values and procedures which cause the wide variation in delivered 
energy of drill rigs is therefore necessary. 

Factors affecting the reproducibility of the Standard Penetration Test include: 
personnel, equipment and procedure. A. Casagrande and L. Casagrande (1968) noticed 
considerable differences in penetration resistance N values obtained by two different drilling 
contractors in sands at the same depth on the same site in Michigan adjacent to Lake 
Michigan. Also from a research study at the University of Colorado (Goble and Ruchti 1981) 
on the performance of safety hammers operated by cathead and rope, it was found that the 
height of the drop as well as the delivered energy varied between different operators and also 
for the same operator. In a comparison study between the cathead and rope system and the 
Mobile Drilling Company’s automatic winch and cable operated driving system, Kovacs et 
al. (1975) performed additional testing to evaluate several parameters that effect the impact 
velocity. These parameters included fall height, hammer type, inclination of the drill stem, 
number of wraps on the cathead, cathead speed, rope age, method of release of the rope, 
sampler type and hammer stroke. Kovacs also presented a useful summary of available 
research on the SPT. 

The effect of in-situ stresses will change the blow count accordingly, but this is the 


in-situ state of the soil at a particular site and the blow count itself is correct. Other factors 
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affecting the blow count and leading to erroneous results are the actual mechanical 
performance of the SPT procedure. These factors are inherent in the present procedure and 
are sometimes difficult to avoid. They were pointed out by several authors and are: 


1. Inadequate cleaning of borehole. 

2. Upward seepage of water into borehole. 

3. The height of drop of hammer. 

4. The friction between hammer and rods during drop. 
5. Weight (type) and length of drill rods. 


6. Monitoring number of blows. 
The most significant of these is probably (3) because this determines the impact velocity of 
the hammer. From the study at the University of Colorado where the velocity of the hammer 
at impact was measured, and the height of drop showed great influence on the velocity. 


Kovacs, et al. (1978) also found similar results and proposed a new testing procedure. 


2.5 Energy Measurements 

In any field or laboratory testing procedure, the ability to reproduce results is 
important. In the case of the SPT, the main source of the unreliability of the blow count is 
the difficulty in producing identical drops to maintain the same impact velocity as well as 
delivered energy. From a study done by Schmertmann and Palacios (1979), it was found that 
the blow count is inversely proportional to the energy absorbed by the soil during a test. In 
order to determine experimental results of this kind, a way of measuring the energy in the 


rods was developed. 
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It appears that the work of Schmertmann et al. (1978) was the first that sought to 

more directly evaluate the variations in ram impact velocity. He made strain measurements 
in the drill rod and calculated delivered energy. In his procedure he assumed that the kinetic 
energy of the driving system is fully transmitted to the drill rod before the reflected stress 


wave arrives back at the top. Thus, the transferred energy is 


a 


E(t!) = = | F(t)°dt (2.1) 


where c is the speed of the stress wave propagation, E is the modulus of elasticity of the drill 
rod, A its cross-sectional area, F is the force measured near the top of the drill rod given as 
a function of time and ¢’ is the time when the first force record passes through zero. The 
analytical basis of Equation 2.1 is presented in Chapter III. 

The primary weakness of Schmertmann’s approach is that for short drill strings the 
reflected wave will arrive at the top before all of the transmitted wave has been input. If 
there is a loose connection in the drill string, an early zero is also possible. 

Kovacs et al. (1978) attempted to measure the velocity of the hammer at the instant 
of impact. This will give a measure of the energy available just before impact occurs. 
Measurements of this type, however, do not include energy losses at impact, and more 
important, do not give a time history record of the event. Therefore, a method from the pile 
dynamics area was used. It is known from pile dynamics that the energy delivered to the rod 


can be expressed as a function of time 
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fu 


E(t) = | F()v@dt (2.2) 


where v is velocity measured as a function of time. These quantities can be measured using 
Strain gages for the force while the velocity can be obtained by integration of the output of 
accelerometers. From this setup a force-dispiacement relationship of the event also would 
be available. This method, however, has been unsuccessful since problems involving the 
acceleration measurement occurred. This was said to be due to high vibration frequencies 
resulting from metal to metal impact (Hauge 1979). 

In an attempt to measure the energy absorbed by the soil during an SPT test, 
Schmertmann and Palacios (1979) used two force transducers, one at the top and one at the 
bottom of the drill string. In this way the incident and reflected energy could be obtained 
from force measurements at the top and the bottom, respectively. The difference between 
the two would then be the energy absorbed by the soil. With this method, a successful way 
of measuring the energy was obtained. Figure 2.7 shows a schematic sequence of the 
hammer-impact wave pulse moving up and down the rods and through the two .oad cells 


mentioned previously. 


2.6 Computation of Bottom Velocity, Displacement and Reaction Force 

From a research program at Case Western Reserve University by Goble, et al. (1975) 
and Teferra (1977), a procedure of finding resistance forces acting at the pile tip during 
driving has been developed. This method makes use of measured force and acceleration at 
the pile top as input quantities. The computational procedure is based on the superposition 


of different boundary conditions. A free pile is defined to have the prescribed measured 
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velocity (or acceleration) at the top and zero forces acting along its length and toe end. The 
actual pile is defined as a pile which has the measured velocity at the top and the real 
resistance forces acting along its length and toe end. Thus, the difference between the 
measured top force of the actual pipe and that of the free pile will be the top force for a pile 
whose top is fixed, and to which the actual resistance forces are applied. This difference has 
been referred to as the Measured Delta Curve. Here, the effect of the actual resistance forces 
on the force at the pile top has been separated from the forces due to the applied velocity. 
In this analysis the skin resistance forces are assumed to be zero, which actually corresponds 
better to the SPT than to pile driving because the hole is augered down to the desired depth 
in the SPT procedure. Therefore, all side friction should be eliminated. 

In order to obtain the Measured Delta Curve an expression must be obtained for the 


free pile solution. This can be obtained from Equation (A.14) in Appendix A. 


F(t) = “00 - wit - *) + v(t - et (2.3) 


Here F(t) is the free pile solution, v,(t) is the measured top velocity, L/c is the time required 
for the wave to travel the length of the pile, and A is the cross-sectional area of the pile. 
Thus, with the aid of Equation 2.3, a curve of the free pile solution is available. The 
difference between this curve and the measured top force is the Measured Delta Curve which 
represents the effect of the tip reaction forces on the pile top force. In the case where no 
skin friction forces are available, as in the SPT procedure, the reaction at the toe can be 
computed from the Measured Delta Curve using Equation A.16 in Appendix A, by letting x, 


equal LZ (i.e., the load acts at the toe of the pile). The following expression is obtained: 


F() = 2R¢ - 2) - 2R¢ - 74) + are - 4 - are - % ~~. (2.4) 
Cc c c Cc 
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where F/t) is the top force effect (Measured Delta Curve) and R(t-L/c) is the toe reaction. 


By rearranging this equation the expression for the toe reaction is obtained. 


L 


reo 2 Arey res 24 
c 2 


peng< 24 ert lS sk (2.5) 
Cc c Cc 


The bottom velocity is computed by superimposing the measured top velocity produced by 
the impacting hammer and the velocity effect produced by the toe reaction forces, assuming 


fixed top and a free toe pile model. 


2 _k,_¢ = _ 3l, , ¢ _ 2b 
v,(f) = 2v¢ =) 7 eae 2(t a wre = 


tL, (2.6) 
se? 


5L ¢ 
+ 2v(t - —-) - —_2R¢ - 
Now, with an expression for the toe velocity as a function of time, the toe 


displacement can be computed by the integration of the toe velocity. 


rf 


dt) = | v,(t)dt (2.7) 


where d,,, is the toe displacement and the time ¢ is measured at the bottom of the pile. 


toe 

With records of the toe displacement, Equation 2.7, and toe reaction force, Equation 
2.5, a means of determining the force-displacement relationship, which yields the spring 
constant of the soil, can be obtained. The method described above was applied to the pile 
driving procedure, but since, as mentioned before, the main assumption for this derivation 
is the assumed zero resistance force along the pile (i.e., only toe resistance), the method 


seems to have the ability to apply well to the SPT procedure. A basis for such an analysis, 


however, is accurate records of the measured force and velocity at the top of the rod. 


CHAPTER III 


WAVE PROPAGATION IN THE SPT 


3.1 One Dimensional Wave Mechanics 

When a rod is struck at the end by a mass a small zone of material is first 
compressed. This compression causes a strain, ¢, and, therefore, a force F = eAE (A is 
cross sectional area and Z is Young’s Modulus of the rod). The force, F, then compresses 
a neighboring particle. However, since material is compressed a motion of the particles is 
also necessary, and a particle velocity, v, is generated in the rod. Whenever a velocity is 
given to a particle of mass m within a time period Ar the particle has to be accelerated 
causing an inertia force mv/At. This inertia force is in balance with the strain force and 
because it takes time to accelerate the particles, the strain will be transferred at a certain 
speed, c, called the wave speed. From the assumption that the wave propagates in the 
longitudinal direction only, (i.e. the transverse effect in a rod, where the ratio of radius to 


length is small, will be insignificant) the governing equation for the wave transmission is, 


from (Goble 1988) 
au 28 G.1) 
Derivation and solution of this equation and a discussion of the boundary conditions are given 


in Appendix A, adapted from (Rausche 1988), (Goble 1988) and (Fischer 1984). 


The speed of wave propagation is only dependent on the material properties of the 
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rod, and not on boundary conditions or the physical behavior of the hammer impact. It is 
equalto /E/p  , where p is the mass density of the material. For the AW steel rods used 


in this project, c was taken as 16,800 ft/sec. 


3.1.1 Boundary Conditions 


In a uniform, elastic rod, where no external forces act, the stress gradient will travel 
through the rod without being changed in magnitude. In Equation (3.1) an infinite rod was 
assumed for the derivation. An important consideration is to deal with what happens when 
a stress wave arrives at the end of a rod. Two cases will be dealt with, namely a wave 
approaching a free and a fixed end. 


a)Free End: 


When a stress wave arrives at an end, the stress gradient will be changed. At a free 
end, the particles will be subjected to higher accelerations since no further material is 
strained in front of the wave. Now, due to higher accelerations, a new stress gradient will 
build up between particles next to the end. A new stress wave will thus be generated 
traveling away from the end, and is called a reflected wave. The particles will move in the 
same direction, while the sign of the stresses will change. Thus, a compression stress wave 
arriving at a free end will be reflected back as a tension wave, having the same absolute 
magnitude as the compression wave, and the same form. Since the stress wave is not 
instantaneous, but has a finite rise and decay time, a portion of this reflected tension wave 
is superimposed back on the still oncoming part of the compression wave. Because the 
velocity in the oncoming and reflected wave has the same direction, the total velocity in the 


immediate region of the tip will double during the time of reflection. 
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b)Fixed End: 


For the case of a wave approaching a fixed end, opposite conditions result. 
Acceleration of the particles at the tip is not possible, and thus the superposition of motions 
at the tip will result in zero velocities. The stress wave will reflect with the same sign, and 
the stress will build up to twice its original magnitude at the point of reflection. A 
compression wave will reflect as a compression wave, and therefore the direction of particle 


motion will be opposite in oncoming and reflecting waves. 


3.1.2 Discontinuity in Bar Properties 

A sudden change in the bar properties will impose a disturbance in the wave 
transmission. Depending on the nature of the change, for example, if an increase or decrease 
in area is encountered by the oncoming wave, part of this wave will be transmitted and part 
will be reflected at the point of change. Two continuity conditions are imposed at the 
juncture, namely, that the force and displacement at both sides of the juncture must at all 
times be equal. From these facts the reflected and the transmitted wave can be calculated in 


terms of force and particle velocity, (see Appendix A). 


= —__— = 3 . 2 
F, 7 “= >a5" (3.2) 
= - 1 = =~ 1 = 3 . 3 
F, 7 F, v, i y, (3.3) 


Where, 
vy, and F, are the velocity and the force in the transmitted wave, 
vy, and F. are the velocity and the force in the reflected wave, 


y, and F, are the velocity and force in the incident wave, 
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a is the ratio of the rod mechanical properties and is equal to Z/Z,; Z, = E,A,/c,, where 
indices 1 and 2 are for oncoming and transmitted section, respectively. 

These relations show that the transmitted wave always has the same sign as the 
incident one. The sign of the force and the velocity in the reflected wave is dependent on 
the material ratio a. For a wave approaching a decrease in area, a < 1, and the sign of the 
reflected stress will be opposite to that of the incident, while the velocity will be in the same 
direction. When a wave approaches an increase in area, ~ > 1, the opposite condition 
results. For a = 1, that is, there is no change in the relative material constants, there is of 
course no reflection and the transmitted wave equals the incident. It is here also interesting 


to note that in the extreme cases when a = 0 or infinity, we have reflection at free and fixed 


end respectively. 


3.2 Proportionality of Force and Velocity 

At impact between hammer and rod a stress gradient is introduced at the end of the 
rod, producing a stress wave that travels down the rod. In a uniform rod this stress gradient 
will cause the same particle velocities independent of the location on the rod. One must here 
distinguish between particle velocity and stress wave velocity. By particle velocity reference 
is made to the motion of a particle point in the rod during the passage of the stress wave. 
This stress wave has a velocity of propagation referred to as c. From wave mechanics one 
can develop a relationship between the stress at a point on the rod and the particle velocity. 
In Figure 3.1 the wave has already traveled to a point which is still at rest. During a time 
interval At the wave travels additional distance AL = c\At and since the material below 


became compressed, the point moved a distance 5. This deformation 6 was caused by a strain 


e over a distance AL and therefore 
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_ 86 
Ee=— 
AL 
But replacing AL by c\At leads to 
e= 
cAt 


Since the point traveled a distance 6 during a time Ar it had a velocity 


o 
y= — or o = vAt 


At 


Substituting for 6 in e we find 


Thus, the strain at a point in the rod material is proportional to the particle velocity of the 


same point. The relation can be expanded to cover stress 


and force 


F = ytA (3.4) 
c 


The proportionality constant EA/c is commonly referred to as the rod impedance because it 
is that force with which a rod opposes a sudden change of velocity by one unit. 

Equation (3.4) is valid for an infinitely long uniform rod or, applied to our case with 
finite rod length. It is valid when a single wave passage is considered. This means it is only 
valid within the first 2L/c time, or the time required for the stress wave to travel down the 


rod, reflect at the bottom and then return to the point under consideration. 


strained 


AL 


unstrained 


Wave Speed 
c 


Figure 3.1: Proportionality of Force and Velocity. 
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3.3 Generation of Waves in the SPT Procedure 

When the hammer impacts on the rod, it creates a compression wave traveling down 
the rod, and at the same time, a stress wave is generated in the hammer. Since the force and 
the velocity are proportional, the shape and the magnitude of the stress waves can be 
determined by looking at the velocity. Fairhurst (1961) describes what ideally happens when 
two similar materials impact. Across the plane of contact two conditions must be fulfilled 
during impact: 

(i) the force in the hammer must equal the force in the rod 
(ii) the absolute spatial velocities of the striking end of the hammer and the struck end 
of the rod must be equal at all times that the two surfaces are in contact. 


From these conditions the particle velocity in the rod and the hammer can be written in terms 


of the impact velocity. 


v, -( “ J» (3.5) 
1+a 
" -( J J» (3.6) 
l+a 
where, v is the impact velocity, 


v, is the particle velocity in the hammer, 
v, is the particle velocity in the rod, 


a is the material ratio. 


3.3.1 Generation of Waves in the Rod 
Consider now a hammer blow in the SPT procedure, using a safety hammer and AW 


rods. The hammer and the driving rod are shown with their dimensions in Figure 3.2. The 
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Figure 3.2: Safety Hammer and Drill Rod, with Dimensions. 
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AW rods come in five foot segments and are connected by threads to make up the total length 
required for the borehole to be tested. The schematic of a blow is shown in Figure 3.3. 

In (a), the hammer is lifted up to a height of 30 inches then dropped onto the top of 
the driving rod, (b). Impact between hammer and rod occurs at point A. A compression 
stress wave is generated in the rod, propagating downward with speed c and two stress waves 
are generated in the hammer emerging from this point. One is going upward and is a 
compression wave, the other traveling downward creating tension in the hammer. The 
distance from the point of impact, A, to the top of hammer, B, is very short and since a free 
end reflection takes place at B, the stresses in section A-B will at most times be canceled due 
to superposition of waves. For convenience and practical purposes, one can thus ignore the 
part A-B and say a tension wave is generated, traveling down the hammer. The wave 
propagation in the system is shown in Figure 3.4. The hammer and the drive rod are of same 
material, which means that the material ratio a@ equals the ratio between the areas, which in 
this case is 0.3. Then from Equation (3.6) the particle velocity in the drive rod y. would be 
0.77v (v is the impact velocity). Just before impact, the particle velocity in the hammer is 
equal to the impact velocity and at impact this velocity will decrease by 0.23v (Eq. 3.5) to 
make the particle velocity in the hammer equal to that in the rod which is 0.77v (Fig. 3.4b). 
The same values would represent the stresses. When the downwards moving tension wave 
reaches the bottom of the hammer, D, a free end reflection takes place. The additional mass 
at the bottom of the hammer is very small compared to the hammer weight and the stress due 
to the inertia forces of this mass is negligible. This tension wave reflects as a compression 
wave of same magnitude which cancels the oncoming tension wave (Fig. 3.4c). At the 
juncture E between the drive rod and the AW-rods, a decrease in area is encountered by the 


oncoming wave. Again the two rods are of the same material, so the material ratio a equals 
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Figure 3.4: Wave Propagation in the Safety Hammer System. 
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the ratio between the areas, which is 0.457, This results in a reflection and transmission of 
waves according to Equations (3.2) and (3.3), which yields 63% of the compression wave 
force transmitted and 37% reflected back up as tension. From Equation (3.2) the transmitted 
particle velocity is 1.373v,, where v,, in this case, is v, which equals to 0.77v (Fig. 3.4d). 

The transmitted wave propagates down the rod and encounters the transducers located 
six inches below juncture E. At the joint between two AW-rods, a new change in area is 
encountered. This change, however, is so small that the disturbance resulting from it have 
minimal effect on the wave transmission. When the downward moving wave reaches the 
bottom of the rod (Fig. 3.3), it is reflected according to the boundary conditions, and travels 
back up the rod. It is again sensed by the transducers, and then it comes to juncture E where 
again a change in area is felt (in this case an enlargement). Therefore a reflection and 
transmission takes place. The transmitted wave proceeds to the top of the driving rod where 
it is reflected back down the rod. At this time an impact is also felt by the hammer. The 
cycle is then repeated. 

The initially downward moving compression wave was transmitted and reflected at 
juncture E. The reflected wave is in tension and its value from Eq. (3.3) is 0.373, where 
y, equals to 0.77; then its value in terms of the impact velocity is 0.287v. This tension wave 
moves up the rod until it is reflected on top. Before this wave reaches the top and at 21,/c 
time (Z, is the hammer length), the compression wave moving up the hammer would reach 
point A where reflection and transmission will take place. The wave transmitted to the drive 
rod is in tension and its value from Equation (3.6) is 0.1480 and the wave reflected down the 
hammer is also in tension and its value is 0.080 (Fig. 3.4e). Shortly after that the wave 
reflected at the juncture E will reach point A and again transmission and reflection take 


place. From Equations (3.5) and (3.6), the reflected and the transmitted waves are both in 
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tension and their values are 0.1850 and 0.102¢, respectively (Fig. 3.4f). The upward and 
the downward moving tension waves superimpose on each other to give a total value of 
0.620. The original downward moving compression wave(0.77¢) will be reduced by 0.620 
and the result is a downward moving compression wave of value 0.150. Again transmission 
and reflection takes place at the juncture E. The resultant transmitted wave is still in 
compression of value (1.373)(0.150)= 0.2060. The downward moving tension wave (. 1850) 
is reflected as a compression wave of value 0.07¢ which moves up the rod until it is 
reflected on top as a compression wave. When it reaches the juncture E it repeats the cycle 
of reflection and transmission as described earlier. 

Since at the top, and at the bottom, one can see irregularities in the hammer 
geometry, many small reflections and transmissions will take place in these regions. Due to 
these facts, a thorough description of the wave propagation in the hammer is very 
complicated and frustrating. Superposition of waves will in some instances cancel their 
effects, end a consideration of waves generated in the hammer longer than the 2LZ/c time of 
the hammer is believed to be unnecessary. Reflections and impacts wil] also cause losses in 
wave intensity, and the overall effect of waves in the hammer will therefore diminish quite 
rapidly. 

The signals sensed by the transducers will at all times be dependent on the 
superposition of waves passing this section. A blow will first be sensed when the wave 
transmitted down through section E passes the transducers. Now, depending on the length 
of rods, a reflection from the bottom or from the top will get to the location of transducers, 
and their intensity will be superimposed on the ones existing. A theoretical estimation of the 
wave sensed by the transducers in a 20 ft long rod is shown in Figure 3.5 along with a plot 


of stress and velocity waves from experimental data. 
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Figure 3.5: Theoretical Wave Form Against Experimental Data. 


The force will oscillate about its static value, with smaller and smaller intensity after 
the impact forces cease, and gradually damp out. For an actual rod in the SPT, the velocity 
will oscillate about the value prescribed by Newton’s second law since it is not fixed on top, 


but moves some distance into the soil. 


CHAPTER IV 


MEASURING AND PROCESSING EQUIPMENT 


4.1 Introduction 


The SPT is a dynamic soil exploration test, and the understanding of the SPT requires 
an understanding of its dynamics. This involves an investigation of the dynamic interaction 
between hammer, drill rod, sampler and soil. The standard penetration number is dependent on 
the energy delivered to the sampler (Schmertmann and Palacios 1979). This energy is determined 
by the velocity of the hammer at impact, and can be determined, if measurements are made on 


the drill rod by Eq.2.2, 
E(t') = f F(t)v(t)dt 


Force and velocity are denoted by F and v, respectively, where the measurements are made above 
the ground surface as a function of time during a hammer blow. The force and velocity describes 
the dynamic motion of a point on the rod. The knowledge of the motion at one point can be 
transferred to another by the aid of one dimensional wave mechanics. The requirements for 
getting a quantitative insight of the dynamics of the SPT are to have accurate records of the force 
and velocity as a function of time at a point on the rod. By measuring these quantities in the rod, 
and not the velocity of the hammer at impact, an accurate measurement of the phenomenon is 
obtained, and it is not necessary to account for energy losses resulting from the hammer impact 


and rope friction, heat loss at impact and so on. 
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4.2 Transducer Description 

From the pile dynamics area the capability had been developed to measure force and 
velocity in the pile during driving. Goble and Ruchti (1981) found that the output from the load 
cell does represent not only the forces in the pile, but also the inertia forces due to the load cell 
weight. Goble et al. (1975) developed two types of force transducers; one type introduces a force 
measuring device between the pile and the hammer, and the second type can more accurately be 
described as a strain transducer. Good accuracy was achieved by mounting strain transducer and 
accelerometer near the pile top. Because of the similarity between pile driving and the SPT 
procedures, it is assumed that the same measuring system can be transferred to the hammer-rod 
system used in the SPT. The measuring system used for piles is described by Goble et al. (1975) 
and the system used on the SPT system will be described here. 

For the force measuring device four 350 Ohm strain gages were mounted directly on a 
five foot long rod about six inches from the top. Gages were mounted at 90 degree intervals to 
minimize local bending and to obtain the average force in the rod. Each gage is a two-element 
gage, one element is parallel and the other is transverse to the axis of the rod. The gages were 
then arranged into a full bridge and directly calibrated using a bridge amplifier. Schematic of 
the gage hook-up and the equivalent bridge are shown in Figure 4.1. The rod was then placed 
in the MTS testing machine for static calibration (Fig. 4.2) and the results shown in Figure 4.3. 
There is a good linear response within the range tested. A strain gage signal conditioner was 
used for the amplification of the signals coming from the gages during the calibration. 

The acceleration of the rod was measured by mounting two accelerometers diametrically 


Opposite on the rod at about the same level as the strain gages. In this way any bending effects 


a) Arrangment of Strain Gages. 


b) Equivalent Bridge. 


Figure 4.1: Strain Gages Layout. 


Output 
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g Device Using the MTS Testing Machine. 


Figure 4.2: Static Calibration of the Force Measurin, 
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Figure 4.3: Static Calibration Results of Force Measuring Device for AW-rod. 
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during a blow could be canceled by taking the average of the two acceleration signals. 
Commercially available Entran accelerometers were used. The type used is EGCS-240D. The 
nominal range of this device is from 5g to 5000g for 15 volts excitation and a frequency range 
of 30% to 50% of the resonant frequency which is 7 kHz. The calibrations of the accelerometers 
were performed by Entran and the results supplied by them, The technical specifications for the 
accelerometer are shown in Table 4.1. The rod was shaped in such a way that the accelerometers 
can be glued to the rod itself on a flat surface. 

With these two electronic devices mounted on the rod a fairly convenient and portable 
way of measuring the dynamic response of the SPT system was available. The rod could be 
taken out in the field during a regular test and mounted on top of the string of drill rods used, 


or it could also be used for experimental testing in the laboratory. 


4.3 Data Processing 

The electronic signals coming from the transducers were amplified and filtered by running 
them through the strain gage amplifier. The excitation voltage was always set to the same voltage 
as that used during the calibration performance and, therefore, no changes in the calibration 
constant is required. The output signals from the amplifier can be recorded by the two data 
acquisition systems used in this project. 

The acquisition systems used are the Computerscope ISC-16 manufactured by R.C. 
Electronics, Inc. and the Pile Driving Analyzer (PDA) manufactured by Pile Dynamics, Inc. The 
Computerscope ISC-16 data acquisition system is a fully integrated hardware and software 
package that consists of a 16-channel A/D (Analog/Digital) board, and external Instrument 
Interface and Scope Driver software. The package is capable of receiving up to 16 channels of 


data input at an aggregate sampling rate of up to 1MHz. The hardware is installed into a 
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Table 4.1: Specifications for Model EGCS-240D Entran Accelerometer, Supplied by Entran. 
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COMPAQI PC with graphics monitor where the signals could be viewed and then saved on a 
floppy disk in digital form. A processing program can now use these records to calculate force, 
velocity and energy in the rod. A flowchart for the program is shown in Figure 4.4. The actual 
program is written in FORTRAN 77 with interactive data input. In order to clarify the 
computation of these values the processing program will be discussed in detail. 

The digitized data on disk is already sorted out in such a way that each channel is written 
in a separate column. Amplifier setting and calibration constants for force and velocity must be 
given first to scale the records. If two accelerometers were used, the average of the two is 
calculated. Now the zero levels of the two signals must be determined. Some fluctuation in the 
Static state is always present in the transducers, so this variation is averaged for both signals and 
subtracted from them. When the zero level is determined, the acceleration can be integrated to 
obtain the velocity. For this process, the trapezoidal method was used. Now using Eq. (2.2) the 
energy can be computed as a function of time. 

An example plot of the results obtained from this program is shown in Figure 4.5. The 
velocity is multiplied by ZA/c to give the velocity in force units. The record shown was taken 
on a 20-foot long rod in the laboratory. As can be seen, the graphical plot gives a better 
description of and insight into the dynamic event of a hammer blow and is a good visual aid in 
the comparison. 

The only disadvantage of the R.C. data acquisition system is the inability to save every 
blow during a regular SPT operation. But for research work in the lab where the hammer is 
lifted and dropped manually the R.C. is very useful. For field testing the PDA is more practical. 

The PDA was designed originally for pile driving analysis, but because of the similarity 
between pile driving and the Standard Penetration Test the PDA was used to take measurements 


during SPT performance. The model used in this project is the GCX. This model has the ability 
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Figure 4.4: Flowchart for Processing Program. 
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Figure 4.5: Plot of Force, Velocity and Energy of a Record Taken in the Laboratory Using the 
R.C. Data Acquisition System. 
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to save every single, second, fifth or tenth blow for up to a total of ninety blows if two channels 
are used or eighty blows otherwise. The records are saved in a digital form. The saved blows 
can be sent to the PC then saved on a floppy disk. The PDA does all the required analysis; it 
scales the force and the acceleration, averages the accelerations (if two channels were used), 
integrates the acceleration to get the velocity and calculates the energy in the rod either by 
integrating Fv, F* or v’ depending on the selected method. For this analysis force and velocity 
calibration, impedance ratio EA/c, and //c should be calculated and set properly. Figure 4.6 
shows a plot of force and velocity for a record taken on a 20-foot long rod recorded on the PDA. 
The smoothness of the curves is due to the low pass filter of the PDA. 

Force and velocity signals can be viewed on an oscilloscope connected to the PDA; also 
a printout of certain values such as maximum energy, maximum force, maximum velocity etc. 
can be obtained. The PDA provides analog output of force, average velocity, average 
acceleration and velocity, and acceleration from each accelerometer. This analog output can be 
recorded on an FM magnetic tape recorder, and later processed and analyzed in a more 
comfortable environment. This has the obvious advantage that records all of the signals, not only 
force and average velocity, of successive blows taken in the field can be automatically stored. 
In this way the dynamic event can be recreated in the laboratory. 

The tape recorder used for this project was a TEAC R71 with a capacity of recording 
seven separate channels, simultaneously. Six channels were used for recording of the dynamic 
event and one channel for fluter compensation; one of the six channels was also used for voice 
notations. Two channels were used for accelerations, two channels for force and velocity, and 
the other two channels for force and energy measured by the SPT Calibrator (discussed later.). 
Any technical information such as the set for each blow, length of rods in hole and site 


description was recorded on the voice channel for ease in the reproduction process. 
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Figure 4.6: Plot of Force and Velocity Record Taken Through the PDA. 
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The Standard Penetration Test Calibrator is a complete system, manufactured by Binary 
Instruments, Inc., to measure force and energy in the drill rods during the SPT. The Model 102 
SPT Calibrator was used in this project. It consists of a 40,000 pound capacity load cell and an 
instrument to compute the energy from the dynamic force and the properties of the drill rods, 


according to the relationship 


E(t) = vi f 2 Ry?dt 


where 

E(t) = energy in the rod. 

A = area of rod cross section. 

c = velocity of the compression wave in the rod. 

E = Young’s modulus for the rod. 

F@ = measured axial force in rod. 

t = time from hammer impact. 
The displayed energy is computed as a percent of the theoretical 4200 inch-pounds available from 
the impact of a 140 pound weight falling 30 inches. The load cell was attached to the drill rod 
at about 2 feet from the point of hammer impact. 

The velocity of the hammer was measured by radar and plotted by the Hammer 
Performance Analyzer (HPA). A picture of the radar and the HPA is shown in Figure 4.7. An 
output from the HPA taken during a field testing is shown in Figure 4.8. The velocity of the 
hammer starts from zero and increases until it gets its maximum value at impact. The peak of 
this plot is the impact velocity. The calibration of the HPA is 2 ft/sec per division. Therefore, 


the impact velocity shown by this plot is 11.9 ft/sec. 
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Figure 4.7: Radar and the Hammer Performance Analyzer. 
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4.4 Dynamic Calibration 

The calibration of the force measuring device was determined under a static load where 
a good linear response was obtained; and the calibration of the accelerometers was given by 
Entran. Having this information it was thought that the rod was ready for testing. However, the 
force and velocity were not proportional due to an error in the calibration constants. The 
question was; which of the signals is incorrect if not both? This question can not be answered 
easily from measurements taken during an SPT blow due to the complicated wave behavior in 
the SPT system. To be able to answer this question the instrumented rod was tested in a simple 
dynamic test. The test and the results will be described and discussed in detail. 

The apparatus is shown in Figure 4.9. The instrumented rod was connected to another 
two five foot rods to make a total of fifteen feet long rod and then it was suspended by four wires 
(two at 8 inches and two at 150 inches from the impacted end). A four foot long AW-rod was 
also suspended by four wires and used as a pendulum hammer. An adjustable guide was used 
to connect the drill rods and the hammer to the wires. The hammer and the rods were aligned 
carefully to ensure uniform contact across the striking surface. The signals were taken through 
the R.C. data acquisition system and the impact velocity was measured by the radar. 

The hammer rod was lifted about one foot then allowed to fall freely and impact the rod 
creating two compression stress waves traveling in the rod and in the hammer simultaneously. 
Since the hammer and the rod are made of same material and have same cross-sectional area, the 
impedance ratio # equals one. Equations (3.5) and (3.6) show that the particle velocities as well 
as the stresses in the rod and the hammer are equal to half the initial values (i.e half the impact 
velocity). When the stress wave reaches the free end of the hammer it reflects back as a tension 
wave with the same magnitude as the incident compression wave which reduces the stresses to 


zero. At the moment the reflected wave comes back to the point of contact with the drill rod, 
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the hammer will be unstressed and a reduction in stress between hammer and rod results. Since 
the waves in the hammer and the rod are equal the stress in the rod becomes zero and thus, a 
rectangular wave of length 2Z,/c is generated in the rod, propagating with the speed c (L, is the 
hammer length). The phenomenon is shown in Figure 4.10. The wave then reflects from the 
free end of the rod and after a time 2L,/c is again sensed by the transducers (L, is the length of 
the drill rod measured from the gages). The reflected stress wave is a tension wave and shortly 
after it is sensed by the transducers it reflects from the impacted end of the rod as a compression 
wave. This compression wave cancels the arriving tension wave and after a time 1/c (twice the 
length from the transducers to the impacted end of the rod divided by c) the strain gage senses 
a zero stress and stays unstressed for 2L,/c (until the entire arriving wave passes), then senses 
the tail of the compression wave. 

A different phenomenon is observed for the velocity wave. The velocity reflects from 
the free end with the same sign and magnitude. Then after 2L,/c the velocity sensed at the 
accelerometer equals the incident velocity or one half the hammer impact velocity. It 
immediately reflects from the impacted end with same sign and magnitude. After 1/c the 
transducer senses double the velocity for a time equals to 2L,/c then senses the tail of the 
reflected wave before it drops to zero. A theoretical form of the stress and velocity waves that 
are sensed by the transducers is shown in Figure 4.11. 

The theoretical impact velocity due to a fall of 12 inches is 8.02 ft/sec (calculated from 
v = V2gh). But for accuracy, the velocity used to calculate the waves in the rod was taken 
from the radar measurements. Figure 4.12 shows the output obtained from the HPA. The 
impact velocity for this blow is about 8.05 ft/sec. The records from the transducers and the 
theoretical form of the stress and velocity waves are shown in Figure 4.13. Velocity1 and 


Velocity2 are the velocities computed from the accelerations taken from the two accelerometers 
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Figure 4,10: Wave Propagation in the Pendulum Test Performance. 
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Figure 4.13: Stress and Velocity Records in Comparison with the Theoretical Form. 
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for the same blow. These results are obtained by using the static calibration constant for the 
force and the calibrations provided by Entran for the velocity. But as can be seen that the 
stresses as well as the velocities are higher than the theoretical waves which are calculated using 
the measured impact velocity. Therefore, new calibration constants are required. The calibration 
for the force was decreased by 3.4% to become 3600 kips/volt instead of 3727 kips/volt. The 
calibrations for the accelerometers were increased by 15% and 27% to become 0.0447 and 
0.0456 mv/g instead of 0.038 and 0.0333 mv/g, respectively. The new calibration constants were 
used to calculate the stresses and the velocities of the same blow and the plot is shown in Figure 
4.14. There is a very good proportionality between force and velocity; also a very good match 


with the theoretical form. 


4.5 Frequency Response 

The metal to metal impact in the SPT produces signals with some very high frequency 
content. In order to reproduce the actual event, the frequency response of the measuring and 
processing system as a whole must be such that it can handle the significant part of these 
frequencies. The sensitivity and frequency response of the measuring system is determined not 
only by the transducers but also by the strain gage amplifier and the data acquisition system. The 
strain gage amplifier used has a low pass filter of 10kHz, while the frequency range of the PDA 
is 1.2 kHz. This means that the signals taken through the PDA will have a maximum frequency 
content of 1,2 kHz, and the PDA will therefore act as a low pass filter. The frequency contents 
of signals taken through the RC data acquisition system depends on the amplifier filter selection, 
and it can be as high as 10 kHz. From (Oran 1988) one finds that, in order to represent the 
actual analog curve exactly with a digital record, it is necessary to take samples of the analog 


signal at a rate greater than twice its highest frequency. The highest sampling frequency of the 
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Figure 4.14: Stress and Velocity Waves Using the New Calibration Constants. 
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RC data acquisition system is 1MHz which is 100 times the highest frequency of the analog signal 
(10 kHz). 

The calibration curves of the accelerometers (See Figures 4.15 and 4.16) show that for 
the calibration factor provided by Entran the useful frequency is 1.2kHz, and after that the 
calibration factor becomes a nonlinear function of the frequency. Providing the calibration curve 
plotted up to 10kHz this function can be determined and used to process signals of 10kHz 
frequency contents. These signals can be handled by using the RC data acquisition system. To 
determine the calibration function the accelerometer is considered to act as a single degree of 
freedom system of mass m, stiffness k, and viscous damping coefficient, c. The natural 
frequency and the damping ratio of the accelerometers have been determined by matching 
theoretical output to the experimental data. Figures 4.17 and 4.18 show a plot of the 
experimental data provided by Entran and the best match of the theoretical output as well as the 
natural frequency and the damping ratio for this match. Derivation and discussion of the 
calibration factor is given in detail in Appendix B. 

Since the calibration factor is a function of frequency, the acceleration, which is a 
function of time, must be transferred to the frequency domain. -For this process the Fast Fourier 
Transform has been used. Now, having both the calibration function and the acceleration in the 
frequency domain, the acceleration can be scaled; then transferred back to time domain using an 
Inverse Fast Fourier Transfer. A processing program similar to the one discussed before is used. 
A flow chart of the program is shown in Figure 4.19. Results obtained from this program are 
shown in Figure 4.20. The plot shows a comparison between the velocity that sas been 
calculated using the calibration factor provided by Entran and the one calculated from this 


program. The corrections that have been made on the calibration factor has a slight effect on the 
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velocity signal and thus, the error in using the constant calibration factor provided by Entran up 


to 10 KHz can be neglected. 
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Figure 4.15: Calibration Curve of the M437 Accelerometer. 
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Figure 4.16: Calibration Curve of the M438 Accelerometer. 
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Figure 4.17: Comparison Between Theoretical and Experimental Data of the M437 Calibration 
curve. 
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Figure 4.19: Flowchart for the Frequency Domain Processing Program. 
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CHAPTER V 


EVALUATION OF RESULTS 


5.1 Laboratory Tests 
With the measuring system described in Chapter IV, dynamic measurements can be made 
on the SPT during a regular test. Before such measurements were taken in the field, the system 
was tested in the laboratory. A drawing of the laboratory set up is shown in Figure 5.1. The 
instrumented rod was connected to three five-foot rods to make a total rod length of twenty feet. 
At the bottom, a six inch long penetrometer was connected to the rod and it rested on a two-inch 
thick piece of plastic. This makes the length below the transducers 20 feet. In addition to the 
top measurements, bottom reaction and displacement were also measured. Strain gages mounted 
on the penetrometer were used to measure the bottom reaction and an LVDT, connected to the 
bottom, was used for displacement measurements. This system is described in greater detail by 
Chen (1990). The signal read by the bottom strain gages is the difference between the reflected 
tension wave and the incident compression wave. In the case of a free rod, the reflected tension 
wave would be of the same shape and magnitude as the incident compression wave and, thus, the 
signal from the bottom transducer would be zero. But in this case where a bottom resistance 
exists, the difference between the incident and the reflected waves is the soil resistance and, 
therefore, the reading from the bottom transducer is the bottom resistance. 
In order to keep the rod standing after a hammer blow, its top was passed through a hole 


made in a steel plate that was held by two guides (Figure 5.1). The drill rod and the safety 
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Figure 5.1: Drawing of the Laboratory Testing Setup. 
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hammer described in Chapter IIT were used in these tests and the measurements were made by 
the RC data acquisition system. The hammer was lifted to the prescribed height and then 
dropped, and measurements of force and acceleration were recorded and processed. Figures 5.2 
and 5.3 show the results from two consecutive blows. The velocity is multiplied by the rod 
impedance, EA/c, in all figures to give the velocity in force units. Both blows are plotted to just 
past the 4L/c time of the rod. Each 2L/c time interval including the beginning of the record is 
marked on the time axis. 

Focusing now on Figure 5.3, the proportionality between the two curves up to the 2L/c 
time is good. A rapid rise in both force and velocity is observed, but the rise time of the velocity 
is greater than that of the force. This indicates that the frequency content in the force signal is 
higher than that in the velocity signal. The average magnitude of the first peak is about the same 
as the theoretical value, 26.4 kips. The length of this peak is the 2L/c time of the drive rod and 
after that time it decreases to about 8 kips and then oscillates about this value to the 2L/c time 
of the rod. At that time, the force rapidly becomes negative due to the reflected tension wave 
coming from the bottom while the velocity shows an increase in magnitude as expected. The 
same comments also apply to Figure 5.2. 

From the energy curve, it can be seen that the energy delivered to the rod reaches it’s 
maximum at the 2L/c time. The magnitude is 298 lb-ft or 85% of the 350 Ib-ft theoretical energy 
calculated from the ram potential energy before dropping. The plot also shows that 70% of the 
maximum delivered energy is transmitted after the 2L/c time of the drive rod. Calculation of the 
energy absorbed by the soil needs additional measurements of force and velocity at the bottom; 
such measurements and calculations are discussed later. 

With these measurements the bottom resistance, velocity and displacement can be 


calculated using the methods described in Section 2.6. This type of calculation was done using 
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Figure 5.2: Laboratory Test Results, Test # 1. 
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the record shown in Figure 5.3. The free rod solution was obtained by substitution of the 
measured velocity in Equation 2.3. It was plotted on the same graph with the actual measured 
force in Figure 5.4. The difference represented by the shaded area in Figure 5.4 is the Delta 
Curve, F(t), described in Section 2.6. Now, by substituting the values of F(t) into Equation 2.5, 
the bottom reaction is obtained and the plot is shown in Figure 5.5 along with the measured one. 
The spikes in the calculated bottom reaction are due to the difference in the frequency contents 
in the force and the velocity signals, but otherwise the two curves follow the same path. After 
calculating the bottom reaction, the bottom velocity can be calculated from Equation 2.6 then 
integrated to obtain the displacement (Eq. 2.7). The plot of the calculated and measured 
displacement along with the bottom velocity are shown in Figure 5.6. The displacement curves 


do not follow the same path, but the final displacement compare very closely. 


5.1.1 Energy Absorbed by Soil 

At any point on the rod, at any time, the force and the velocity in the rod have the value 
of the instantaneous sum of all incident and reflected waves reaching that point at that time. The 
top transducers provide a record of the incident compression wave, but only until its cutoff by 
the reflected tension wave shortly before 2L/c. At the bottom, the force is measured by using 
the strain gages mounted on the penetrometer and the velocity is measured by mounting an 
accelerometer on the rod at one foot from the bottom. Shortly after L/c time the bottom force 
transducer provides a record of the subtraction of the reflected tension wave from the incident 
compression wave while the accelerometer provides the addition of the reflected and the incident 
waves. By subtracting the bottom force record from that of the top, one obtains the force in the 
reflected tension wave and by subtracting the top velocity record from that of the bottom, the 


velocity in the reflected wave would result. Now with the force and velocity records from the 
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Figure 5.4: Comparison Between the Actual Rod and the Free Rod Solution. 
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Figure 5.5: Comparison Between Calculated and Measured Bottom Reaction. 
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top transducers one can use Equation 2.2 to calculate the incident compression wave energy 
entering the rod. And after calculating the force and the velocity in the reflected wave, one can 
use the same equation to get the energy in the reflected tension wave. By subtracting these two 
energies, the energy absorbed by the soil is obtained. 

Figures 5.7 through 5.11 show the plots for an example of applying the analysis described 
above. Figure 5.7 shows the force and the velocity records obtained from the top transducers 
along with the energy in the incident compression wave. Figure 5.8 shows the superposed force 
records obtained from the top and bottom transducers and Figure 5.9 shows the same plot for the 
velocity records. The plot of the bottom records is shifted back by Lic time so that a better 
comparison can be made. The force and the velocity in the reflected wave, obtained by the 
subtraction discussed before, in addition to the reflected energy are shown in Figure 5.10. The 
energy delivered to the soil is obtained by subtracting the energy in Figure 5.10 from that of 
Figure 5.7 and the results are shown in Figure 5.11. In this case, the energy delivered to the 


rods is 85% of the 350 lb-ft theoretical energy, Z,, and the one absorbed by the soil is 42% of 


E,. 


5.2 Field Tests 

After good laboratory measurements were obtained, the system was tested in the field 
during actual SPT operations. Measurements were taken at sites in the Denver area at different 
depths, with different soil conditions and two different hammers. The CME Automatic and the 
Pilcon hammers were used in these tests. The load cell for the Binary Instruments measuring 
system was mounted on top of the instrumented rod which was connected to the drill rod. A 
picture of the data acquisition systems used in the field and described in Chapter IV is shown in 


Figure 5.12. Examples of results from these measurements are shown in Figures 5.13 through 
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Figure 5.8: Superposed Force Records. 
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Figure 5.10: Force, Velocity and Energy Obtained from the Reflected Wave. 
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Figure 5.11: Calculation of the Energy absorbed by the Soil. 
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5.16. The first thing to notice here is the difference in the shape of the waves produced by each 
of the hammers in the field and the laboratory tests, and also the magnitude and length of the first 
peak. These differences are due to the difference in the geometric shapes of the hammers and 
impact block assemblies. 

Figures 5.13 and 5.14 show a plot of the force, velocity and energy for two blows taken 
through the RC during the CME hammer test. The total length below gages for these blows was 
77.5 feet, and the bottom was a soft rock. The SPT blow count (N-value) was ninety blows per 
foot. The effect of the sampler on the wave propagation can be easily noticed by looking at the 
behavior of the force and velocity at the 2L/c time. The sampler has larger area than that of the 
drill rod and thus at the juncture between them a reflection and transmission takes place. Since 
an enlargement in the area is encountered, part of the incident compression wave will reflect back 
as compression, but the reflected velocity wave will have an opposite sign. These reflected 
waves propagate upwards for a 2L,/c time (L, is the sampler length) before any reflection from 
the bottom arrives. This explanation can be seen in the two figures where, at 2L/c time, the 
force increases sharply (indicates compression) while the velocity goes in the other direction. 
The length of these spikes is 2L,/c. Another thing to notice here is that the reflected wave is in 
compression and that is due to the hard material (soft rock) at the bottom. The impact velocity 
was recorded for every blow, and for the blows shown in Figures 5.13 and 5.14 was 12.2 and 
12.1 ft/sec, respectively. Maximum energy calculated from Equation (2.2) (graph) and the 
energy obtained from the Binary Instruments measuring system are shown on each of the plots. 

Figure 5.15 also shows the force, velocity and energy as a function of time for another 
blow taken through the RC during the CME performance tests. The length below gages was 47.5 
feet, and the soil was sandy clay. The SPT blow count was four blows per foot. At the first 


peak, notice that the velocity is smaller than the force and then at about 5 milliseconds the 
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Figure 5.13: Field Test # 10, Blow No. 50, Length Below Gages is 77.5 ft, Taken Through the 
RC During the CME Automatic Hammer Performance, June 1990. 
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Figure 5.14: Field. Test # 10, Blow No. 128, Length Below Gages is 77.5 ft, Taken Through 
the RC During the CME Automatic Hammer Performance, June 1990. 
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Figure 5.15: Field Test # 8, Blow No. 5, Length Below Gages is 47.5 ft, Taken Through the 
RC During the CME Automatic Hammer Performance, June 1990. 
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Figure 5.16: Field Test #8, Blow No. 9, Length Below Gages is 82.5 ft, Taken Through the 
PDA During the Pilcon Hammer Performance, August 1990. 
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velocity goes above the force. This lack in proportionality is believed to be due to using one 
accelerometer only. That was because the connection of the lead wire to the plug of one of the 
accelerometers had been cut off. The effect of using one accelerometer is presented in Figure 
5.17 where for the same blow shown in Figure 5.12 the velocities from each accelerometer were 
plotted separately along with the force. In both plots the proportionality is not as good as when 
both accelerometers were used. 

In Figure 5.15, a different behavior of the force and the velocity can be observed, after 
the 2L/c time, when driving in soft material. The reflected force wave is in tension while the 
velocity is reflected back with same sign as the incident one. 

Figure 5.16 shows a plot of a record taken, at a different site, through the PDA during 
the Pilcon hammer performance test. The length below gages was 82.5 feet, and the soil 
consisted of soft rock. The N-value was one hundred twenty blows per foot. At about 9.5 
milliseconds a sharp drop in the force and a rise in the velocity was appeared. This was found 
to be due to a lose connection at that depth. This lose connection delayed the wave propagation 


by 1.78 milliseconds and made the 2L/c time 11.6 milliseconds instead of 9.82 milliseconds. 


§.2.1 Stress Peak Evaluation 


As can be seen from the previous results, the value of the stresses generated in the rod 
from a hammer blow is dependent on the hammer’s geometric shape. The shape and the 
magnitude of a stress wave generated in the rod by the safety hammer was discussed in Chapter 
II. The maximum stress peak was 1.060,, where o; is the stress obtained from multiplying the 
impact velocity by E/c of the rod. This value is 6% higher than that if the hammer were a rigid 
mass and struck on the drill rod itself. If the anvil was removed and the hammer was dropped 


directly on the drill rod, the maximum stress would be 0.80;. The maximum stress value 
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obtained while using the pilcon hammer was the same as that of the safety hammer, (see Figure 
5.16) but different stresses were obtained during the CME hammer performance. This 
phenomena gave a reason for looking at the wave propagation when the CME hammer is used. 

The CME automatic hammer is a single acting hammer which consists of a 140 lb ram 
that impacts on the top of an anvil connected to the drill rod, The dimensions of the ram and the 
anvil are shown in a drawing of the system in Figure 5.18. The ram and the anvil are of same 
material which means that their impedance ratio equals the ratio between their areas which in this 
case is 0.25. Just before impact the ram is moving as a rigid body with a velocity equal to the 
impact velocity, v. At impact, a part of this velocity is transferred to the anvil and a stress wave 
is generated, propagating downwards in the anvil and upwards in the ram. From Equation 3.6, 
the velocity transferred to the anvil is 0.8v and the one reflected to the ram is 0.2v. The same 
values would represent the stresses. The wave propagation in the system is shown in Figure 
5.19. At the juncture between the anvil and the drill rod, transmission and reflection take place. 
The anvil and the drill rod are also of the same material and thus, the impedance ratio is equal 
to their area ratio which is 0.224. From Equation 3.2 the transmitted velocity is 1.307v and the 
stress is 1.307v. Therefore, this value is 30.7% larger than that if the drill rod was struck by 
a rigid mass. The reflected wave is in tension and moves up the anvil until it is reflected on top. 
The hammer is still in contact with the rod, a transmission and reflection takes place, and a 
reflected tension wave moves down the anvil and repeats the cycle of reflections and 
transmissions as described earlier. At the same time, the wave traveling up in the hammer would 
reflect at the free end and propagate down until another impact is felt by the anvil and repeats 
the cycle. 

If the anvil were removed and the ram were dropped directly on the drill rod, the stresses 


generated from such a blow would be smaller than that if the anvil were used. The impedance 
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Figure 5.18: Dimensions of the CME Automatic Hammer System. 
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Figure 5.19: Wave Propagation in the CME Hammer System. 
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ratio between the hammer and the drill rod is the ratio between their areas which is 0.056. By 
substituting for this value in Equation 3.6, the stresses in the rod would be 0.9476;. 

It is of interest to note that the stress value of 1.3070, is about 30 ksi for the velocities 


measured with the radar. This compares very well with the peak stress shown in Figure 5.17. 


5.2.2 Bottom Analysis 

As mentioned in Section 2.6, one can, from top measurements, calculate bottom 
quantities such as resistance, velocity and displacement as a function of time. These quantities 
give a better understanding of the soil behavior under a hammer blow. To demonstrate this, 
bottom analysis was made for two blows taken in different soil conditions. The first blow is the 
one shown in Figure 5.14 where the bottom consists of soft rock and the second was shown in 
Figure 5.15 where the soil was a soft, sandy clay. Figures 5.20 and 5.21 show the soil resistance 
calculated for each of the blows. In the sandy clay material where the blow count was very low, 
the soil resistance vanishes very rapidly and the sampler moves very easily, this movement is 
represented by the velocity and the displacement in Figure 5.22. The final displacement is 2.8 
inches and that gives a blow count of four blows per foot which agrees with the N-value counted 
in the field. Different behavior is observed in the soft rock strata. The velocity vanishes quickly 
and gives small displacement (Fig. 5.22) while the bottom resistance is large. The displacement 
here is 0.2 inches which gives a blow count of 60 blows per foot whereas the blow count from 
the field was 90 blows per foot. Also, plots of the bottom reaction versus displacement gives a 
better explanation of the event and those are shown in Figures 5.24 and 5.25. 

The other important thing that can be calculated from these bottom quantities is the 
energy delivered to the soil and this is obtained from the multiplication of the reaction by the 


displacement. Figure 5.26 shows the delivered energy as a function of time for the record shown 
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Figure 5.20: Bottom Resistance for the Blow Taken in the Sandy Clay Soil. 
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Figure 5.21: Bottom Resistance for the Blow Taken in the Soft Rock Strata. 
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Figure 5.23: Bottom Velocity and Displacement for the Blow in the Soft Rock. 
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Figure 5.24: Bottom Reaction Versus Displacement for the Record Taken in the Sandy Clay 
Soil. 
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Figure 5.25: Bottom Reaction Versus Displacement of the Record Taken in the Soft Rock. 
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Figure 5.26: Energy Delivered to the Soil for the Record Taken in the Soft Rock. 
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energy delivered to the soil and this is obtained from the multiplication of the reaction by the 
displacement. Figure 5.26 shows the delivered energy as a function of time for the record shown 
in Figure 5.14. The maximum energy is about 175 lb-ft which is 55% of the energy transferred 
to the drill rods and 50% of the theoretical E,,. 

This analysis describes the sampler behavior under a hammer blow in two different soil 
conditions and also provides a force-displacement relationship from which soil properties can be 
obtained. No attempt has been made to study the engineering properties of the soil or its 


correlation to the SPT blow count. 


CHAPTER VI 


CONCLUSIONS AND DISCUSSION 


The experiments and testing described here have given a better insight to the measurement 
of the dynamics of the SPT. With the successful measurements of force and acceleration at the 
top, the door is opened for further analysis which may lead to a more accurate interpretation 
of the SPT in general, and also this might yield valuable information for pile driving analysis 
and capacity. 

The main problem that has, in the past, faced the dynamic measurements on the SPT is 
the frequency response of the measuring system. In this project, the frequency response of the 
measuring system has been shown to be satisfactory. Both signals, the force and the acceleration, 
were filtered to have a frequency contents not more than 10 KHz. The accelerometers response 
for these frequency contents was found to be adequate and the sampling frequency of the data 
acquisition system is high enough to satisfy the needs. The calibration factor of the force 
measuring device obtained from static tests and that of the accelerometers provided by Entran 
were found to be incorrect and thus the system was recalibrated under dynamic load using the 
pendulum test. This test also gave a better understanding of the response of the measuring 
system. 

From studying of the wave propagation in different systems used in the SPT procedure, 
it was found that the stresses generated in the rod from a hammer blow depends on the system 


that is used in this procedure. Stresses generated by using the safety hammer were lower than 
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that generated by the CME automatic hammer, but the difference in the systems seems to have 
no effect on the energy delivered to the rod. Also, changes in the system itself were shown to 
have an effect on the stress level in the rod. Removing the drive rod from the safety hammer 
system or the anvil from the CME system decreased the stresses in the rod, but no further study 
has been made to find if this fact has any effect on the blow count. 

Analysis and calculations of bottom quantities from top measurements for laboratory and 
field tests are presented in Chapter V. This analysis is a demonstration for what one can obtain 


from the top force and acceleration measurements and it can be the key for further studies on 


the SPT of this aspect. 
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APPENDIX A 


ONE DIMENSIONAL STRESS WAVE MECHANICS 


IN A ROD UNDER IMPACT 


A.1 Derivation of the Wave Equation 


In order to study the characteristics of the stress wave, the displacement due to the strain 
associated with the passage of the stress wave will be considered. Let the displacement of a 
particle in a uniform, elastic rod at distance x from a fixed coordinate origin at time f to be 
u(x,t). After a time dt the point at the distance x will move a distance u to a new position X and 
the neighboring point at x+dx will displace by u+du to a position of distance X+ dX (Fig. A.1a). 


The change in length of the portion dx is : 


dX - dx = du = “dy 
ax 


and therefore, 


_ AX - dx _ ou 


£ 
ax ox 


Now, considering the free body shown in Fig. A.1b and using Newton’s second law 


(F = ma) the following relation can be obtained: 
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Figure A.1: Derivation of Wave Equation. 
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oo _ eu (A.1) 
ox or 
But, 
o = E.e and e = ou 
ox 
Therefore, 
Ox ox ax? 


Substitute (A.2) into (A.1), and the governing differential equation, the so-called wave equation, 
is obtained 
Fu _ ru (A.3) 
ar? ax? 
Where c” = E/p. E and p are Young’s modulus and mass density of the rod, respectively. The 


general solution of Equation (A.3) is 


u(x,t) = f(x + ct) + g(x - ct) (A.4) 


which is easily checked by back substitution. The solution u(@v,4) = flx+ct) represents a wave 


traveling in the negative x-direction and u(x,t) = g(x-cf) represents a wave traveling in the 


positive x-direction. 


A.2 Wave Propagation 


The displacement of a point at distance x, at time ¢t, is given by 


u(x,,t,) = fa, + ct,) + gy, 7 ct,) 
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After a time ¢ = t, + Ar the wave propagating in the positive x-direction travels a distance Ax 


while the wave propagating in the negative direction travels a distance -Ax and therefore 


u(x,t) = f(x, - Ax + clt, + At) + g(x, + Ax - cit, + Ad) = (A.5) 
Now, if c is defined as the speed of the wave propagation, then Ax equals cAr. Substitute for Ax 


= cAt in Equation (A.5) and simplify 


u(x, = f(x, + ct.) + g@, - ct,) 


or, 


u(x,t) = u(x,,t,) 
This shows that after time Ar the stress wave, which is represented by the displacement, 
propagates a distance Ax = cAft with a velocity c keeping its shape and magnitude. Notice that 


the earlier definition of c shows that it only depends on the material properties of the rod. 


A.3 Relation between Force and Velocity 


Differentiating Equation (A.4) with respect to time leads to an expression for the velocity 


v(x,t) of a particle 


af(n) _ ate) (A.6) 


- 9 - 
w(x,t) = ag tt) e| an au 


where y = x + ct and p = x- ct. The above equation can be written as 


v(x,t) = vx + ct) + V(x - ct) 


thus, v can also be described by two waves, v, and v,, which travel in opposite directions. 
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Now, suppose that the shape of the function f and g and, therefore, u(x,t) is known at a 


certain time ¢. Then the strain e is obtained by differentiating u(x,t) with respect to x 


ea) = = 2 [fetch + ge - od) 


or, 

e(xt) = An) , ost) _ wy, ~y) (A.7) 
mH op c 8 

The advantage of using derivatives with respect to y or yw becomes apparent since these 

derivatives express the behavior of the derivative with respect to both x and t. Using the rod 


cross-sectional area A and taking compressive forces positive, the force F(x,t) in the wave can 


be calculated. 


F(x,t) = a [ye + ct) - v(x - ct)] (A.8) 


Thus, there exists a simple relation between particle velocity and force in a stress wave; the force 
is proportional to the velocity by a factor AE/c. From Equation (A.8) it can be seen that the force 
will be compressive when the velocities of particles and wave propagation have the same direction 


and it will be a tension force otherwise. 


A.4 Boundary Conditions 


So far, only the homogeneous differential equation has been considered and nothing has 
been said about boundary conditions or external forces. In a uniform, elastic rod, where no 
external forces act, the stress gradient will travel through the rod without being changed in 
magnitude. In the above derivations an infinite rod was assumed. An important consideration 
is to deal with what happens if there are any external forces acting along the rod or when a stress 


wave arrives at the end of arod. The one dimensional wave equation is linear and superposition 
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is valid, so that a case of complicated boundary conditions can be split into several basic types 
of easily solvable problems. If the boundary conditions were nonlinear, then a result can be 
obtained by assuming piecewise linear boundary conditions and superimposing their effects. A 
wave reaching the end of the rod might encounter either prescribed force or displacement 
conditions. This problem can be split into the case where the wave travels in a rod with 
homogeneous (i.e. zero force or displacement) boundary conditions plus the case where no wave 
is present with non-zero end forces or displacements. The case where external forces are acting 
along the rod may also be treated separately and then superimposed onto the homogeneous 
solution. Prescribed displacements have to be considered at the end of the rod only. Therefore, 
four basic conditions must be treated: 

¢ Wave approaching a free end 
¢ Wave approaching a fixed end 
e Prescribed force acting at a point along the rod 


e Prescribed displacement at the end of the rod. 


A.4.1 Wave Approaching a Free End 

Consider Figure A.2a where on the left a stress wave o,(x-ct) is shown approaching the 
free end (x=L) of the rod. The condition of zero stress must be maintained at all times at that 
end. This condition may be satisfied by a second stress wave o,(x+ct) propagating towards the 
left, which, when superposed on the incident wave, cancels the end-section stresses. Expressing 


this concept mathematically 
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Figure A.2: Reflection of Displacement and Stress Waves at Free End. 
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from which 
CA 


Oe (L = ct) 


of, 
ral - et) = - 
Hence it is evident that the slope du/dx of the left propagating wave must be the negative of the 
slope of the forward-propagating wave as each part of the waves passes the end of the rod. The 
displacement waves shown in Fig. A.2b demonstrate this condition, and the corresponding stress 
waves show clearly how the stresses at the tip are canceled. 

Although the concept of a left-moving wave coming from beyond the end of the rod 
makes it easier to visualize the mechanism by which the boundary condition is satisfied, it should 
be understood that this wave actually is created at the end of the bar as the forward-propagating 
wave reaches that point. In other words, the incident wave is reflected at the free end; the 


reflected wave has the same deflection as the incident wave, but the stresses are reversed. 


A.4.2, Wave Approaching a Fixed End 


To consider now the case where the right end at the bar is fixed rather than free, the 


condition of zero displacement must be satisfied. Therefore, the displacement at x = L can be 


expressed mathematically as 


un, =~9 =f, (L - ct) +f, + et) 


from which the reflected wave may be expressed in terms of the incident wave as 


AL - ct) = f,@ + ct) 
Thus the displacement waves in this case are seen to have opposite signs, and by analogy with 
the preceding discussion it can be inferred that the incident and reflected stress wave have the 
same sign, as shown in Fig. A.3. Hence, in satisfying the required zero-displacement condition, 


the reflected wave produces a doubling of stress at the fixed end of the rod. 
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Figure A.3: Reflection of Displacement and Stress Waves at Fixed End. 
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A.4,3 Prescribed Force Acting at a Point Along the Rod 


If a force is applied at x = x’ along the rod, then the continuity condition requires 


u(x',,t) = u',,t) 
where x’, and x’, are on the left and the right side of the loaded cross-leftline, respectively 


(Figure A.4). This can also be written in terms of velocity, 


UE) _ UC pe!) (A.10) 
at at 


Also the equilibrium condition has to be satisfied, i.e., 


AE [See _ air) = F (x',) (A.11) 
ot ot 


where F,,(x’,t) is the applied force at x = x’ and time ¢. Now recalling Equations A.6 and A.7, 


Equations A.10 and A.11 can be satisfied by choosing 


of(n) _ Osh) ~ipy 
| an o | au .. 5 FA's t)/EA 


and 
al - eval - 
on =1'p Op a, 
Substituting these two conditions in Equation A.10 leads to 


c [jae eA ‘ =¢ (0 + SFaes0)/ 2A 


thus, satisfying the continuity condition. 
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Figure A.4: Impact Force Applied at an Intermediate Point along the rod. 
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Similarly it proves that the equilibrium condition is satisfied. Figure A.4 shows that both 

waves carry a force of one-half of the applied force while the velocities are the same in both 
waves having the same direction as the applied force and magnitude, 1/2FA(x’,)c/EA. 

The case where the force is acting at the end of the rod can be deduced since one of the 

half-waves is immediately reflected and superimposed on the other. Therefore, a wave will travel 


away from the end with particle velocity, 


u(x,t) _ op ,-b-% 
ot A c 


and strain 


A.4.4 Prescribed Displacements at the End of the Rod 

The case of a prescribed displacement condition at a point along the rod will not be 
discussed since this is equivalent to a rod of shorter length with prescribed end displacement. 
A displacement prescribed at the end of a rod is equivalent to a prescribed velocity, v,(Z,f). 
From the force-velocity proportionality relation this can be considered a force condition where 


the force has to be chosen as F,(Z,f) = v,(L,)EA/c. Thus the previous case can be applied. 


A.5 Superposition of Waves 


In this leftline certain special cases will be treated where the findings of the previous 
sections are applied. These special cases will be needed in treating problems where external 


forces act along the pile together with prescribed end conditions. 


108 
A.5.1 Free Rod Under Known Velocity at the Top 
Suppose a velocity, v,¢), which is zero for t < 0, is imposed on the top of a rod of 
length, L. It is desired to know the top force necessary to maintain this velocity when no other 
forces are acting along the rod. 


As long as the wave created by v,(4), has not yet reached the free end (i.e., for t < L/c) 


a rod particle at a point has velocity, 
*; 
v(t) = v(t - ~ 


At time t = L/c the wave reaches the free end of the rod and the free end reflection case 
applies. Thus, a reflection wave having velocities of the same sign and forces of opposite signs 
will travel back up the rod. The velocity at a station x = x, becomes 


v(t) = ¥, ( - a4 +y, i. - *| (A.12) 


c 


fort < (2L-xJ/c. 

A new reflection wave will be generated when this up-traveling wave reaches the top 
since here the velocity is prescribed. The case of a fixed end reflection describes this fixed end 
situation. This second reflection requires a force; therefore, the proportionality between applied 


velocity and top forces will no longer hold. The force at the top of the free rod may be denoted 


by Ff); then fort < 2L/c, 


F, (t)4y,(z) 
c 
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and for 2L/c < t < 4L/c 


F, (t) = “ty - vce - 24) (A.13) 
c c 


At a time t = 3Z/c the wave will again be reflected at the bottom end. (It had been reflected at 
the top at time 2L/c). At time t = 4Z/c the wave has to be reflected again at the top but this time 


it has the opposite stress due to the previous reflection at the free end. In general, for 2Lr/c < 


ts 2L/ (r+), 
F(t) = 4’ ly +2 Ye - jay (A.14) 
c jl 


where r = 0,1,2,..., refers to the time interval considered. Equation A.18 gives the exact 
solution for a given top velocity and no reaction forces. If this solution is subtracted from the 
measured force, the Measured Delta curve is obtained as defined in Chapter II. 

Equation A.12 gives an expression for the velocity at some point in the rod for times ¢ 


< (2L + x)/c. The equation can be extended for time (Q27 + x)/c St < QL + 1)+x)4, 


r = 0,1,2,..., 
y(t) = ECD! p [i - 7%) -y, i. - gen) | (A.15) 
j-0 c c 


A.5.2 External Force Acting on a Rod with Fixed Top and Free Bottom 
The force denoted by R,(¢), is assumed to act upwards. If it is acting at x = x, then the 
results of Section A.4.3 apply. Thus, two waves are introduced traveling in opposite directions 


and having particle velocities (at x = x,, and time 4), 


v,(f)=- 5 R(t) 
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The forces in the waves are 


1 

F(t) = +—— 
rf 2R(t) 

i.e., a compression in the upwards and tension in the downwards traveling wave. The reaction 

force F,,,(¢) on the fixed top of the rod due to this velocity will be a compression force of twice 


the magnitude of the force in the wave. The upwards traveling wave will arrive at the top at a 


time x,/c after it is applied. Hence 
*; 
Fyp(t) = R(t - 4) 


which is valid for a time as long as it takes the initially downwards traveling wave to reach the 
top after reflection at the bottom. This reflection at the free bottom end causes the initially 
downwards traveling wave to change the sign in force and therefore, is also a compression wave 
after reflection. Reaching the top at a time (2L-x,/c after it was generated by the resistance 
force, this wave also will produce twice the force at the top which it was propagating. The next 
wave to arrive at the top will be the initially upwards traveling wave, its sign will be converted 
so that a tension wave arrives. 
This way the reaction at the fixed rod top due to R,(f) can be calculated for 


0<t< 4LA, 


4L - x, 


x, 2L - x; 2L + x, 
F 7 (4) = R(t - > + Rt - —;? - Rt - —| - R(t - ) 


For later times it can be observed that all the waves which arrived at a time ¢, at the top have in 
the meantime changed their sign twice, thus, arriving again with the same sign at a time ¢t, + 


4L/c. This result can be expected since a pile of length L has a lowest natural frequency of c/AL, 
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inferring that in the absence of external forces a harmonic behavior can be observed. Therefore, 


for times 4r7L/c < t < (7 + 1)4L/c, 


-F _* _ AL _ 2be% AL 
Fip(2) [Ac 2 = 7) + R(t - j42)| 


pip - cet _ AL 4 Ga 4k (A. 16) 
Rc Fo Re a) 


and r = 0,1,2,..., indicates the time interval considered. 

The problem of obtaining the particle velocities at a point x - x, when the load, R,(/), is 
acting at x = x, is somewhat more complicated. It can be split in two parts, first considering the 
wave which is initially moving upwards and then the wave which is initially moving downwards. 
Both of these waves have a particle velocity -(c/2EA)R{?). In order to facilitate the derivation 
it is further assumed that x, > x; Then the velocity v,,,(¢) at x = x, due to the upwards 
traveling wave obtains its first contribution after a time (x, + x,/c, i.e. the time necessary for 
the wave to reach the top and upon reflection the station x = x,. The wave will again be 
reflected at the bottom end with no sign change in velocities and reach x = x, a second time with 


the same sign in the velocities. Thus, when observing the wave’s action at x = x, for a time f¢ 


< 4L/c, 


2L + x, ~ hy 


v,(t) = £ | (¢- 7%) + RG - 
huup = ia |* c R, c 


2L + x, + 


x: 
*) - R(t - 


4L - x, - 
meen) (A.17) 


similarly the influence on the velocity at x = x, of the initially downwards traveling wave can 


be determined. 
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- x, 2b - x, - 
Vado?) = se |- Rit - BE) RG eee) 
. ee 2-5-*) ia 


Observing that the waves arrive with the same sign of force and velocity after every 4Z/c, then 
superimposing the results from Equation A.17 and A.18, the velocity v,(t) due to R,(t) becomes, 


for 47L/c St Ss (r + IAL, 


c x, +X; 40 2L+X,-%, AL } 
tf) =-—-— t- - j— t- -—_———  —- }— 
v,() 2EADy" Ln - J—) + Rf : J ) 

2b. + Xe +X 4L 
-R(t - fa) (bs ie ayae) 
2b =X, =X 
“R(t - At ~ jy - Roe - —— 44 - fy 


21x, ~-%, AL ~~ x AL 
+R RC 


where again r = 0,1,2,..., indicates the time interval considered. 


A.5.3 Free Rod with a Known Force at the Top 

In this case a force F(t) is prescribed at the top and no forces act along the rod. It is 
desired to obtain an expression for the velocity at the rod top v,,(f), due to this force. The 
reflections of waves at both ends of the pile have force but no displacement restrictions. Thus, 
if the applied wave induces compressive strains, then all the reflection waves will have tensile 
velocities. Reflections will always add to the top velocity on every wave arrival at the top. 


From proportionality and due to the discussed reflections the result can be readily obtained. 
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E 4L 
Viop (t) = al F(t) + 2F,(¢ - 2) + 2F,(t - oe | 
or using the notion the above equation can be rewritten to yield for 27L/c s t S (+ 1)2L/c, 


p(t) vn = [A (t) 2 nc - 4) 


jal 


A.5.4 External Force Acting on a Rod with Free Top and Bottom 
The external force may again be denoted by R,(¢) and act at x = x;. The reflections will 
be of such a nature that the waves arriving at the top always add to the top velocity v,,(¢) if R,(@) 


does not change sign. Hence for 0 < t < 2L/c, 


Ripe) SR 
: Cc c 


eee 
EA 


and for 27L/c < t S (7 + 1)2L/c, 


2L - x, 
Vigp (t) = pd Rae - a at - 522) + Rt - are ae 


EA} 


A.6 Discontinuity in Bar Properties 

The wave reflection which takes place at the fixed or free end of a uniform bar may be 
considered as special cases of general reflection and transmission phenomena occurring at any 
discontinuity in the bar properties. The conditions of equilibrium and compatibility which must 
be satisfied at all points along the bar require that additional reflected and refracted waves be 


generated at the juncture between bars of different properties in response to the action of any 
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incident wave. 
Consider, for example, the juncture between Bars 1 and 2 shown in Fig. A.5. The 
properties of the bars on each side of the juncture are characterized by their mass per unit length 


m and axial stiffness EA. Also, the wave propagation velocity on each side is given by 
£ 
p 


A forward propagating wave, u,, which arrives at the juncture in Bar 1 is reflected as u, which 
travels in the negative direction in Bar 1, and at the same time creates a refracted wave u, which 
propagates forward in Bar 2. 

Continuity and equilibrium conditions are imposed at the juncture: 
Displacement: u; = uy, u, + Uu, = U, 
Force: F, = F,, Fi+F, =F, 
Where the fact that both incident and reflected waves act in Bar 1 has been indicated. Because 


this continuity condition must always be satisfied 


ot ot 


Ou, , Ou, _ OU, (A.19) 
ot 


but the incident wave can be expressed in the form 


u, =f - ct) = f(n) 


Now, the derivatives of u, can be expressed as 
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Figure A.5 Wave Reflection and Refraction at Discontinuity. 
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from which it is evident that the time and position derivatives are related by the velocity of 


propagation 
Ste. _g Me (A.20) 
ot Ox 

similarly 
Oy te, (A.21) 
ot Ox 
Ou, coe (A.22) 
Ot ox 


where the positive sign in Eq. A.21 is due to negative direction of the reflected wave. 


Substituting Eqs. (A.20), (A.21) and (A.22) into Eq. A.19 yields 


Ou fe) 
eye 4g Me = Wg, he (A.23) 
ox ox ax 
but the strains 


@ 
—— =e etc. 
ox 


can be expressed in terms of forces acting in the bars 
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G F 
e,=—=—* ete 
E AE, 
Thus Eq. (1.23) can be expressed as 
c c C 
-_— a * —F, b- __F. ¢ 
AE, A,E, A,E, 


or more simply 


where 
_ c,A,E, 
CAE, 
or 
Z. 
a= 
Z, 


-1 

F.=F {=& 
b (SS) 
F,= F, 22) 
a+il 


These relations show that the transmitted wave always has the same sign as the incident one. The 
sign of the force in the reflected wave is dependent on the material ratio a. For a wave 
approaching a decrease in area, a < 1, the sign of the reflected stress will be opposite to that 


of the incident wave. For a wave approaching an increase in area, « > 1, the opposite condition 
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will result. For a = 1, that is, there is no change in relative material constant, there is of 
course no reflection and the transmitted wave equals the incident. It is here also interesting to 


note that in the extreme cases when a = 0 or infinity, we have reflection at free and fixed end 


respectively. 


APPENDIX B 


ACCELEROMETER CALIBRATION FACTOR 


B.1 Calibration Factor 

The EGCS-240D accelerometer consists of four semiconductor strain gages bonded to a 
simple cantilever beam which is end loaded with a mass (See Figure B.1). The strain gages are 
connected to make a fully active Wheatstone Bridge. Damping is achieved by using a viscous 
fluid medium. Under acceleration, the force on the cantilever is created by the acceleration effect 
on the mass (F=ma). The accelerated mass creates a force which in turn provides a bending 
moment to the beam. This moment creates a strain which results in a bridge unbalance. With 
an applied voltage, this unbalance produces a voltage deviation at the bridge output, which is 
proportional to the acceleration vector. 

In order to study the effect of the frequency on the accelerometer calibration factor, the 
accelerometer was modeled as a damped single degree of freedom system of mass m, stiffness 
k, and viscous damping coefficient c. The model is shown in Figure B.2. If x() is the 
displacement of the mass m and y(‘ the displacement of the vibrating body, the equation of 


motion of m is 


mi + c(x - y) + k(x - y) =0 (B.1) 


Letting the relative displacement of the mass m and the case attached to the vibrating body be 


Zz = (x - y) then, 
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Figure B.1: Entran’s EGCS-240D Accelerometer Cutaway. 


Figure B.2: Accelerometer Model. 
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z=(-y) 


and 


Now, by substituting for (x - y), (¢ - y) and # in Eq.(B.1) the following equation of motion is 


obtained 


mZ + y) + cz + kz =0 


or 


mZ + cz + ke = -my (B.2) 


The motion of the vibrating body is considered to be periodic, then y can be expressed 


in terms of Fourier series as 


and therefore, 


substitute for y in Eq.(B.2) 
mz + cz + kz = my w; Y, eit (B.3) 
jel 
Equation (B.3) represents the equation of motion of the mass m with respect to the case 


of the accelerometer. The steady-state solution of this equation is 


j=l 
from which, 
iw,t 
= = 
z iy’ w, Ze iw,z 
j=l 
» _ = 2 iwe _ 2 
z=l yy; Zz = “Wz 
7) 
and 


> ett - 


jal 


Ns 


Substitute for these values in Eq. (B.3) to get 


2, 2Y, 
“mw, Z + icw,z + kz = mw; —Z 
Zz 
i 
divide by mz 
2,,¢€ k 2¥; 
-w, + i—w, + — = wr (B.4) 
i m J m fj Z, 
But, 
xk. w, and <= 2ew, 
m m 


where, w, is the natural frequency 


& is the damping ratio 
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Now, Eq.(B.4) can be rewritten as 


Y, 
2»; 2 j 
“w, + 2ibwiw, + Wo = WS 


WN 


divide by we and let r = w, /w, 


Y, 
_ + DiFp = pe? 
(1 - r?) + 2iér = r? Z 


f] 
or 


Zi r 


¥, (1 - r?) + 2iér 


multiply the numerator and the denominator by {(1 - 7°) - 2itr} and arrange we get 


Zz  Pad-r) _, 2Er3 
a — {—— + 
Y, (1 - 7?) + (2ér)’ (1 - r?)? + (2Er?? 


or 
Z, 
+ =A+iB 
Y, 
where, 
A= rd - r?) 
(1 - r?)? + (2Er)? 
B= 2Er? 


( ~ ? + (28n? 
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Z, and Y, are both complex numbers, then the ratio of their magnitude is 


i = |= = VA? + B? = ee 
yo 4 (1 - 9°) + (28r)? 
Therefore, 
_ 2y _ 2y 
Z, = rs RY, _ Rwy; (B.5) 
VQ - 72)? + (Er? Ww, 
where, 


ee 
(L - 7°)? + (Er)? 
wi, is the amplitude of the acceleration of the vibrating body and w2 is constant because w, is 
a property of the system; therefore the relative motion z is proportional to the acceleration y if 
the magnification factor R is constant for all ranges of operation. 
A periodic vibration generally has a number of harmonic components, each of which 
gives a corresponding value of r. Amplitude distortion occurs if the magnification factor R 
changes with the harmonic components. It is evident then that the parameters involved are the 
frequency ratio r and the damping factor £. Figure B.3 shows a plot of R for various values of 
damping ratio £. The diagram shows that the useful frequency range of the undamped 
accelerometer is somewhat limited. However, with £ = 0.7 the useful frequency range is 0 < 


wiw, <= 0.2 with a maximum error less than 0.01 percent. 
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Figure B.3: Plot of the Magnification Factor R for Various Values of the Damping Ration é. 
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B.2 Correction 
When the frequency of the vibrations to be measured is higher than that of the useful 
frequency of the accelerometer, the magnification factor R becomes a nonlinear function in the 
frequency ratio r. Thus, the relative motion z becomes non-linearly proportional to the 
acceleration. 
To get the actual acceleration of the vibrating body a different procedure must be used. 
The measured displacement z is a complex function in ¢ and by using Fast Fourier Transform 
(FFT), it can be transferred to another complex function in w, For each value of w, the 
magnification factor R can be evaluated and then the acceleration can be computed from Eq.(2.5) 


Z,(w;) 


j=1 to N 
R(w,) 


¥,(w,) = 


where, N is the number of sampling points. 
Now, the acceleration can be transferred to the time domain using inverse Fast Fourier 


Transform. The results are in complex form, but the interest is in the magnitude which is the 


real values. 


